Curs 5 – Semnale stationare si nestationare Testul unit root Seriile de timp stationare intuitiv inseamna medie si deviatie standard constante in timp In aplicatii insa intalnim de obicei marimi nemarginite ca domeniu de variatie Le numim serii de timp / semnale nestationare Aceste semnale nu pot fi aproximate decat eventual cu ajutorul unui trend liniar De aceea ele trebuie diferentiate de cate ori este nevoie pentru a obtine serii stationare, a caror variatie si predictibilitate caracterizeaza sursa de semnal Referitor la predictibilitate, identificarea si caracterizarea comportamentului unei surse de semnal se reduce la determinarea unui spatiu de probabilitate si a naturii acestuia Prezentam in acest curs notiunile de semnal stationar si nestationar alaturi de testul Dickey-Fuller de tip unit root pentru testarea naturii stationare sau nestationare a unui semnal 1 Stationaritate si nestationaritate Aceste doua notiuni au in mintea multora o semnificatie vaga si de cele mai multe ori sunt folosite impropriu In acest curs vom da definitii exacte asupra acestor concepte si vom prezenta testul unit root (Dicke !"uller# de determinare a stationaritatii unui semnal $ sa incepem insa cu ceva care de o%icei sperie studentii Definitie: $ σ-algebra peste o multime & este o colectie Σ de su%multimi ale lui & care este inc'isa su% operatiile de complementare si reuniune numara%ila( )# Σ ≠ φ *# Daca A ⊂ Σ atunci si &+A ⊂ Σ ,# Daca A)- A*- A,- ⊂ Σ este o secventa numara%ila de su%multimi din Σ- atunci A) ∪ A* ∪ A, ∪ - ⊂ Σ /xemplu( &01a-%-c-d2 Atunci Σ 0 1φ- 1a %3- 1c d2- 1a % c-d2 2 este o σ!alge%ra 4ultimea partilor unei multimi este tot o σ!alge%ra 5erec'ea (&- Σ# pentru care Σ este o σ!alge%ra a lui & se numeste spatiu masurabil- adica se poate defini o masura peste Σ care ataseaza un numar real pozitiv numit masura oricarei multimi din σ!alge%ra Σ De ce σ-algebre 6 5ro%a%ilitatile de evolutie ale unei surse de semnal se definesc de o%icei pe secvente de stari sau pe multimi integra%ile 7e%esgue in spatiul starilor 8otiunea de σalgebra sta la %aza notiunii de camp de pro%a%ilitate care descrie universul de evolutie al sursei Cel mai 9rau: caz care se poate intampla atunci cand analizam sursa unui semnal- este ca ea sa ai%a un comportament pro%a%ilistic- adica sa nu stim ce va emite la momentul urmator de timp decat dupa o anumita distri%utie de pro%a%ilitate Aceasta distri%utie de pro%a%ilitate se poate defini nu doar peste starea curenta a sistemului In cel mai 9rau: caz sursa de semnal are si 9memorie: de capacitate 8- iar iesirea ei depinde de intrarile cat si de iesirile inregistrate anterior pana la maxim 8!) pasi in trecut Distri%utia iesirii unei astfel de sursa de semnal este deci definita peste un spatiu & x ; de dimensiune 8 sau iesirile trecute σalgebra este cea care generalizeaza notiunea de dependenta a iesirii de 9ceva:( de numere (o singura intrare reala#- de vectori de intrare- de secvente de intrare 8!dimensionale (sistemul cu memorie din primul exemplu# si in general peste orice multime 7e%esgue integra%ila Spatiul de probabilitate ataseaza pro%a%ilitati iesirii unei surse de semnal in functie de orice element al unei σ-algebre: numarsecventa- cu%- cur%a- in general orice multime integra%ila 7e%esgue Asadar de acum incolo cand discutam de identificarea comportamentului unei surse de semnal si nu stim nimic despre aceasta- vom lucra la nivel intuitiv cu σ-algebra multimilor integrabile Lebesgue peste Rn, unde n 0 numarul intrarilor sursei @S A ! serie de timp nestationara- si seria variatiilor zilnice (dail returns# ;eturns?n305rice?n3!5rice?n!)3 – serie stationara 8$TA( /volutia seriilor de timp nestationare nu poate fi estimata dincolo de o regresie liniara decat pe portiuni si decat pe termen scurt- si doar daca distri%utia de pro%a%ilitate a sursei generatoare de semnal are o variatie lenta in timp 5rezentam in continuare testul unit root Dicke !"uller pentru investigarea naturii stationare sau nestationare a unei serii de timp estul unit root Dic!e" #uller In spiritul o%servatiei anterioare- cele mai uzuale aproximari ale seriilor nestationare sunt modelele liniare auto!regresive (A;# @n model auto!regresiv de ordin p- notat A;(p# este( ()# unde ϕ)- ϕ*-ϕ,- ϕp sunt parametrii modelului si ε este zgomot al% (spectrul de frecvente uniform- DT"T 0 constanta# Consideram asadar un model A;()# simplu pentru sursa generatoare de semnal( ?n3 0 a ?n!)3< εn (*# unde a este o constanta si ε un proces sto'astic de medie = Introducem operatorul 9diferenta inapoi:( ∇ ?n3 0 ?n3! ?n!)3 Avem( ∇ ?n30 (a!)# ?n!)3 < εn ∇ ?n3 0 δ ?n!)3 <εn (,# - sau (B# Dic!e"-#uller test )( Daca δC= atunci spunem ca seria de timp nu are unit root si este stationara Altfel seria este nestationara 4ai exact( - Daca δ 0 = atunci spunem ca se satisface ipote$a null si ( ?n3 0 ?n!)3< εn si seria se numeste random walk nestationara - Daca δD= atunci seria de timp este doar nestationara Dic!e"-#uller test * ! @nit root cu drift( 5resupunem ca avem un drift constant in variatia lui ?n3 (adaugam a= in ecuatia (,##( ∇ ?n3 0 a= < δ ?n!)3 <εn (5# Testul * testeaza existenta unui unit root cu drift - Daca δ0= spunem ca seria de timp este random Ealk cu drift a= - Daca δC= seria este stationara cu drift - Daca δD= seria este nestationara cu drift Dic!e" #uller test %: Testul , verifica ipoteza ca sursa sa emita un semnal cu o componenta de tip trend determinist Astfel ipoteza este ( ?n3 0 a ?n!)3 < % n < a= < εn (F# Componenta %Gn este componenta de trend pentru semnale discrete- sau %t – componenta trend pentru semnale continue Aplicand operatorul 9diferenta inapoi: avem( ∇ ?n3 0 a= < δ ?n!)3 <%Gn < εn Dupa estimarea parametrilor- daca %≠=( - Daca δC= atunci seria de timp se numeste trend-stationarasi pentru a o%tine o serie de timp stationara tre%uie scazut termenul %Gn- operatie numita detrending - Daca δ0= atunci seria de timp se numeste random &al! cu trend- este nestationara si are unit root 5rintr!un eventual detrending nu o%tinem decat un random Ealk- deci operatia nu intereseaza decat in cazul δC= - Daca δD= atunci seria de timp se numeste nestationara cu trend De ce este util testul Dicke !"uller6 Analiza "ourier pe semnale nestationare nu este utila 5roprietati interesante ale semnalelor se determina pe semnale stationare Aplicam ∇ de cate ori este nevoie cf testului D" si eliminam driftul sau trendul pana cand semnalul devine stationar fara trend si fara drift 'nexa )# @n camp de probabilitate este un triplet (&- Σ -5#- unde & e o multime- Σ este o σ!alge%ra peste & si 5 este o functie care asociaza oricarei multimi A din Σ un numar real pozitiv mai mare sau egal cu = numit pro%a%ilitate Cea mai simpla si intuitiva Hustificare a faptului ca un camp de pro%a%ilitate este definit peste o σ!alge%ra se regaseste in urmatoarele argumente( a# Daca A⊂ Σ atunci &+A ⊂Σ (cerinta ca σ!alge%ra sa fie inc'isa su% operatia de complementare#- pentru ca mereu sa ai%a loc relatia 5(A# < 5(&+A# 0 ) %# Daca A)- A*- A,- ⊂ Σ sunt multimi disHuncte- atunci( 5(A) ∪ A* ∪ A, ∪ - # 0 5(A)# < 5(A*#< () )um se construieste * + Daca STAT/ e domeniul semnalelor care definesc starea sursei la un moment de timp si $@T5@T este domeniul de variatie a iesirii sursei de semnal- atunci & 0 STAT/ x $@T5@T 